Spin-Charge Recombination and Superconductivity in the t — J Model 
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A spin-charge binding mechanism is analyzed using the Schwinger boson representation of the 2-D 
t-J model. In the normal state spinons are paired into singlets which evolove out of the antiferro- 
magnetic phase at half-filling. The singlets induce, via the hopping term, a pairing of holons and 
bind with the latter. The d x 2_ y 2 symmetry of the order parameter and the shape of the T c vs doping 
curve, observed in cuprate superconductors, emerge naturally as a consequence of the symmetry of 
the spinon pairs, thus connecting superconductivity with quantum spin fluctuations in the insulator. 
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The normal state of cuprate superconductors is essen- 
tially two-dimensional, has a large Fermi surface but ap- 
parently hole-like charge carriers and shows non-Fermi 
liquid charateristics. Anderson has argued that this is 
due to spin-charge separation in a doped quantum an- 
tiferromagnet In this picture, charge is carried by 
spinless holes of concentraion S moving about in a quan- 
tum liquid of spin-singlets. The electron itself is not a 
coherent excitations. The singlets evolve out of the un- 
doped phase. Their presence is consistent with the ap- 
pearance of a spin-gap in the normal state and may 
have something to do with superconductivity, as many 
now believe. However, the electron excitations in the su- 
perconducting state seem to be coherent, appearing as 
sharp peaks in the angle-resolved photemission spectrum 
||. Yet, as recent experiments show [0, the "quasipar- 
ticle" peaks do not exist above T c . Hence, if the normal 
state has spin-charge separation, superconductivity must 
be accompanied by spin-charge recombination. 

Based on the 2-D t-J model, a spin-charge binding 
mechanism for superconductivity was originally proposed 
in 1992, in which spin singlets induce a pairing of holons 
and bind with the latter in order to recover the kinetic en- 
ergy lost in forming the magnetic moments || . Because 
of the pivotal role played by the singlets, it is essential 
that the theory be consistent with the behavior at 5 = 0. 
Such a description was developed using the Schwinger bo- 
son representation. This choice is natural since a mean- 
field approximation gives a good description of Neel or- 
dering and the singlets at half- filling (6 = 0) ||. The 
singlets are described in terms of a paired spinon con- 
densate with p-wave symmetry, which evolves into the 
doped phase where the holons form a Fermi liquid. 

Here we focus on the superconducting transition and 
show that the nature of the superconducting state is in- 
timately connected with that of the insulator through 
the singlets. We will show that the superconducting or- 
der parameter is constrained to have the observed d x 2_ y 2 
symmetry as a consequence of the (p-wave) symmetry of 
the paired spinons. The observed shape of the T c vs dop- 
ing curve (showing a maximum) occurs naturally in our 
theory, and at small 8, T c oc e -1 / 5 , also due to symmetry 
and statistics. The theory predicts the opening of a gap 
at T c which has interesting observable consequences in 



the underdoped region. 

Our approach is based on the following. (1) Spin- 
charge recombination occurs directly through the hop- 
ping term, and in both the normal and superconducting 
states. The energetic reason is as follows. At half-filling 
one has only spins (local moments). A "spin" is local- 
ized charge, and costs kinetic energy of order t, which 
the system tries to recover upon doping. Since kinetic 
energy gain is clearly the largest for free electrons, the 
best way to recover it is to bind spin and charge into 
electrons and restore the Fermi liquid. But to move co- 
herently, electrons must go around each other. This will 
happen for large dimensions (e.g., d > 3), but clearly 
not in d = 1. The situation is intermediate in d = 2 
(2) Recombination is a collective phenomenon, and 
thus is missing at the mean-field (i.e., Hartree-Fock) level. 
It is taken into account by a random-phase approxima- 
tion (RPA) in which the electron appears as a collective 
excitation. There are some misconceptions about the 
Schwinger-boson representation because the Fermi sur- 
face consists of pockets in the mean-field approximation 
0. By including RPA we have shown that, spin-charge 
recombination in the (d = 2) normal state is strong 
enough for electrons to appear deep in the Fermi sea, 
where the kinetic energy gain is large. This is sufficient 
to restore the large Fermi surface and destroy long-range 
magnetic order. However, close to the Fermi surface elec- 
trons remain incoherent (5). (3) Then the singlets give 
rise to superconductivity by binding with holon pairs to 
recover additional kinetic energy. This process is taken 
into account by anomalous RPA, as shown below. A 
similar mechanism, but with an on-site interaction, has 
been recently studied in the SU-2 version of the slave- 
boson theory |1. However, the slave-boson approach is 
rather different because of different statistics, the addi- 
tional symmetry and a very different normal state since 
recombination processes described above are not consid- 
ered. 

The t — J model is described by the Hamiltonian 



H 



(1) 



subject to the constraint that no two electrons can oc- 
cupy the the same site. Here the sum is over nearest 
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neighbors, and for cuprates, t/J ~ 3 — 4. The projected 
electron operator c\ a is decomposed as c ia = b\ a hi, 

where b\ a creates a bosonic spinon of spin a and h% 
destroys a fermionic holon. Then the charge and spin 
densities can be written as rij = h\hi, S + = b\*bii 

and Sf — ^[b^b^ — ftjibij,], provided that the constraint 

h\hi + b\ a bi a = 1 is satisfied. We will treat the con- 
straint on the average. However, the representation has a 
remarkable property. Even when constraints are treated 
on the average, the condition that electrons (and Cooper 
pairs) avoid each other is always satisfied since 



AAi = 



(2) 



because holons are fermions (hihi — 0). Therefore, a 
crucial part of the large-?/ physics is automatically taken 
into account in our theory. 

We will treat the exchange term at the Hartree-Fock 
level. But RPA contribution will be included in treat- 
ing the hopping term. This term describes spinon-holon 
interaction: 



e ( k ) c L c k<x = ^X] e ( k ) 6 Uki,>ki4 2 6 k+k2,<n ( 3 ) 



where the potential e(k) = —2t(c,osk x -l-cosfe^) is simply 
the free-electron hopping energy. As shown earlier |7j] , in 
the mean-field approximation holons form a cosine band. 
The mean-field spinon Hamiltonian is given by 

k,cr k 

where Cfe(k) is associated with hopping. The second term 
describes singlet physics through spinon pairing with an 
anomalous self energy 



(k) = 2JA{smk x +sinfc y ) 



(4) 



where A is the magnitude of the singlet pairing operator 
\\bi-\bji — biibji]. The symmetry of <fi(k) is determined 
at half-filling and remains unchanged in the doped phase. 
As far as magnetic properties are concerned, spinon pair- 
ing is analogous to superconducting pairing, and leads to 
a spinon pseudogap and a suppression of the uniform 
magnetic susceptibility below the pairing temperature 
kTsG ~ J A. Diagonalization leads to a spinon spectrum 
which, near its minimum at ko, is approximately given 
by Wfc(k) = [A 2 + c 2 (k - ko) 2 ] 1 / 2 , where c s is related to 
the spin-wave velocity, and A s = u> m i n > 0. Bose con- 
densation at ko gives rise to magnetic order (LRO) at 
Q = 2ko. At half filling, Q = (n, it), but is incommensu- 
rate away for 5 > 0. Generally A s = in the condensed 
phase, but becomes finite once LRO is destroyed. 

As shown in ref. (5), LRO is destroyed by spin-charge 
recombination in the normal state, which is taken into 
account by RPA. The mean-field electron Green's func- 
tion G c o is the convolution (bubble) of spinon and holon 



Green's functions: G c0 (p) = (A^/3) -1 Y^ q G b (q)G h {p + q), 
where /3 = 1/kT and p = (k, u), lo is the Matsubara 
frequency. In the absence of LRO, G c o is incoherent [||. 
Summing the bubbles, we obtain the the RPA expression 
for normal-state Green's function 



G cn (p) = 



G c0 (p) 



l-e(k)G c0 (py 



(5) 



This is characterized by a large Fermi surface and has 
poles for k deep in the Fermi sea, but none close to the 
Fermi surface |3|. 

Thus the low-energy physics of the normal state is de- 
scribed by a Fermi liquid of holons which have with their 
own Fermi surface centered at (7r,7r), and a condensate 
of paired spinons. It is useful to note that since G c q is 
a convolution, a gap in either the spinon or the holon 
spectrum will induce a gap in G c q so that 7mG c o = 
below some frequency. This will in turn lead to gap in 
G cn since from Eq. (5) ImG cn — 0, unless G cn has poles 
in the gap region. In our case, there are no poles in the 
low energy region. 

Superconductivity involves another type of spin-charge 
binding which requires summing anomalous bubbles. We 
use functional integral method involving three fields: b pa , 
hp and an auxilliary fermionic field a pa which is intro- 
duced through a Hubbard-Stratonovitch transformation 
on the hopping term, as described in ref. (5). The hopping 
part of the action is then given by: 



Ahopp — ^ ' £ 1 s ) a pa- a pcr + [ 



^pcj^pG 4~ C.C. 



(6) 



where the electron field is the composite object c pcr = 
(NP)^ 1 / 2 h* q b p + q . Therefore the term in the bracket 
describes the interaction. The remaining part of the 
action is quadratic. The a-field is conjugate to the o 
field and their Green's functions are simply related. Let 
G a {p) =< ap<j a p<j > an d F a {p) =< a_ p |a p | > be the nor- 
mal and anomalous Green's function, respectively. Simi- 
larly, for the other fields. Then electron Green's function 
are found from the identities: G a (p) = e(k) + e 2 (k)G c (p), 
F a (p) = e 2 (k)F c (p). Hence the introduction of a-fields 
allows us to treat spinons, holons and electrons in a sym- 
metric fashion. Physically G a is the renormalized hop- 
ping potential. Thus, the bare function is: G a o = e(k). 

Next we compute the self energies self-consistently to 
leading order, which corresponds to Hartree-Fock and 
RPA. We use the subscript n and A for normal and 
anomalous functions, respectively. Since the three fields 
appear symmetrically the self-energy of one field is simply 
the bubble of the other two. Thus the normal self-energy 
S a „ is the bubble G c q. This gives 



Gan(j>) 



e(k) 



l-e(k)G o0 (p) 



(7) 



which is the screened hopping potential, and leads to the 
RPA Green's function (Eq. 5). 
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Now, since spinons are paired Fb(p) —< b- P ib p ^ > is 
nonzero in the normal state. Then, if we integrate out 
the spinons we obtain an interaction of the form B 

J2 a h a -p'lK-p h -q+P' Fb ^ + C - C - ( 8 ) 

Clearly such an interaction encourages pairing. We 
can compute the relevant anomalous self-energies ("gap 
functions") by allowing the holon pairing amplitude 
Fhip) =< h-php > to be nonzero (which ensures that 
F a ^ 0) and doing a pairwise decomposition of Eq. (8). 
The result is 

^aA( P ) = ^J2 F h(P^(P-P^ (9) 

Pi 

X hA (p) = 1 ^J2F;(p 1 )F b (p l -p), (10) 
Pi 

E 6A (p) = 0(k) +4g£f.(Pl)^(p-Pl) f (11) 

pi 

where, for completeness, we have included the additional 
contribution to the Y,bA due to holon and electron pairing. 
However 0(k), which is nonzero in the normal state, is 
the driving term, i.e., we seek solutions such that all three 
Y>a vanish if 4>(k) = 0. The Green's functions are to be 
calculated self-consistently. They have the same general 
form 

F(p) = H A (p)/[(G- 1 (p)G- 1 (-p) ± E* A (p)i: A (p)], (12) 

where the minus sign in the denominator is for Ft,. The 
inclusion of the anomalous bubble Y, aA in the Green's 
functions is equivalent to summing these bubbles to all 
orders, thus taking into account spin-charge recombina- 
tion in the superconducting state. 

Note that, if Y, A {-p) = ±E A (p), then Eq. (12) 
yields F(—p) = ±F(p). (Here — p = (— k, — ui)). Since 
the holons are single-component fermions, Fh(—p) = 
—Fh{p). Now, the driving term </>k is odd. Let us there- 
fore assume the full T,bA is odd. Then it follows from Eqs. 
(9-10) that T, aA (p) is even, T,hA(p) is odd. Using these 
in Eq (11), we find that "EbA is odd, thus completing the 
self-consistency loop. 

The holon pairing problem is analogous to the BCS 
problem with spinons playing the role of phonons. As a 
first approximation we specialize to the low-frequency or 
or weak-coupling (i.e., BCS) limit. This will enable us 
to extract some important results that are qualitatively 
correct. First we ignore the second term in Eq.(ll) so 
that YjbA ~ </>(k). Next we put uo — in self-energies and 
define the gap functions A 0) /,(k) = T, ai h A (k,ui — 0). The 
effective action for the holons can be written as 

A e h ff = X> - e h (k))h;h p + ME) p + h - v h p )l 



where e/j (k) is the effective holon energy. This is formally 
identical to the BCS problem. The holons are gapped 
with energy E h (k) = [e|(k) + Al(k)] 1 / 2 . The difference 
here is that A^ (roughly ~ kT c ) is determined by holon 
pairing as well as spin-charge binding. 

At low frequencies, Fb can be approximated by its 
mean- field form: F(,(k, uS) — 4>(k)/[ui 2 +w^(k)], but with 
a renormalized spinon energy w^k. Since spinons are 
pseudogapped, we put w = and replace Wb(k) by a con- 
stant Wo which will play the role of a frequency cut-off. 
This yields the coupled gap equations 

A h (k) = J-j ]T 0(q - k)A a (q).g Q (q), (13) 
q 

Aa(k) = iV^ E ^ ~ q)A.(q)^(q), (W) 
q 

where f3g a , h (k) = ^ ^ h (k, w )/A 0ih (k). 

Using (f>(k) — 2JA(sink x + sink y ) and the fact that 
that A Q (k) is even and A^(k) is odd, we find from Eqs. 
(13-14) that the gap functions must have the form 

A a (k) = ai cos k x + a-i cos k v , (15) 

Ah(k) = hi sin k x + h% sinfcj,, (16) 

where, cij and hi are constants. On general symmetry 
grounds the magnitude of A a ,h should be invariant un- 
der k x «-> k y . This implies that A a (k) = ai(cosfc x + 
e l9a cosfcj,), and A^(k) = hi{sa\k x + e l6h sinfc a ), where 
9 a and Oh are unknown phases. Finally, using the fact 
that in our theory < c^c^ >= is always satisfied, we 
have 

y^ai(cosk x + e l9a cosk y )g a (k) = 0. (17) 

k 

Now, g a (k x ,k y ) = g a {k y ,k x ), which implies 6 a = n. 
Then we have 

A Q (k) = -^-(cosfca; - cosfcj), (18) 

where A a o = 2ai . Noting that the pairing amplitude for c 
and a-fields are proportional, the order parameter for the 
physical electron pair has (i-wave symmetry, solely as a 
consequence of statistics and symmetry of the anomalous 
spinon self-energy </>(k). There is thus a deep connection 
between quantum spin fluctuations in the antiferromag- 
net, the spin-gap phenomenon and the superconducting 
state. 

Using Eq.(18) in Eq.(13) we find that Oh = n, i.e., 

Aft(k) = ^-(smk x - sink y ), (19) 

where A^ = 2h\. Hence, the holon gap has p-wave 
symmetry. The gap equations for A a o and A^o are then 
given by 
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A a0 = 



JAA 



Nu; 2 



,i0 E 

k 



(sinfcz -s,m.k v fg h {\s), (20) 



A 



hO 



JAA aQ 



E 

k 



(cos k x - cos k y ) 2 g a (k). (21) 



The fact that the symmetry is determined in the weak- 
coupling approximation does not limit its generality since 
by continuity the strong-coupling function would have 
the same symmetry. 

We can also obtain an approximate formula for T c as 
follows. At T c we can put A a = in the denominator 
of the Green's functions and do the frequency sums [|l0| . 
Combining Eqs. (20) and (21) we get the condition for 
T c : 1 = J 2 A 2 



l2 It%/u)'o, where 



d 2 k (sin k x 



4ir 2 



2e h (k) 2kT c 



d2k ( L. 



cosfc y ) 2 e 2 cr (k) 



(22) 



(23) 



Here e scr (k) = G an (k, lu = 0) is the screened hopping en- 
ergy (Eq. 7) . The quantity tn is a renormalized hopping 
parameter and is independent of T c . The integral / de- 
scribes the physics of holon pairing and, and apart from 
the symmetry factors, has the standard (BCS) form. For 
small 6 (i.e., in the region of interest) holons form a small 
Fermi circle of momentum kph — ('ind) 1 ^ 2 . Then we can 
approximate sin k x by k x , etc, and use the quadratic form 
for the holon energy: e^(k) = t] x (k 2 — k 2 F ), where th is 
the renormailized hopping parameter for holons. Doing 
the integral / we obtain 



kT c = 



_ 2u) X e< 1/x 

IT 



(24) 



where 7 is the Euler constant, u>i is a frequency cutoff, 
and 



t 2 T 2 
tRj -A 2 5 



Uath 



(25) 



is the dimensionless coupling constant that determines 
T c . We can take uj\ = mm(wo, tFh), where tph = iirthS 
is the holon Fermi energy. A more elaborate choice does 
not change the result qualitatively. Based on previous 
analysis 0,0], we expect the parameters tu, th and ujq to 
vary slowly with S compared with the factor A 2 S in Eq. 
(25). The spinon pairing amplitude A is largest at half- 
filling, and decreases with increasing 5 (already at the 
mean- field |j) level. The crucial factor of 5 arises from 
the symmetry factor (sink x —smk y ) 2 in Eq. (22) which, 
at the holon Fermi surface, scales as ~ (k x — k y ) 2 
And. Hence, T c has a maximum, and in particular, 



k 2 - 



at small 6, which is a consequence of symmetry and 
statistics. Thus the characteristic shape observed in 
cuprates occurs naturally in our theory. Furthermore, 
the exponential dependence on A implies that T c is ap- 
preciable only in a small finite region along the 8 axis, 
as observed. These results are qualitatively correct, for 
quantitative accuracy one needs to solve the strong cou- 
pling equations (9-11). 

The holons acquire a gap below T c . As discuused 
earlier (see below Eq. 5), this leads to a gap in the elec- 
tron spectrum. Note that this gap is independent of any 
pseudogap that may exist in the normal state. In other 
words, if the electron already has a gap above T c , the 
total electron gap will show a sharp increase below T c as 
temperature decreases. Apparently such an increase has 
been seen j|. 

To conclude, we have considered a spin-charge binding 
mechanism which differs from traditional approaches in 
that the electron itself is a collective object, and that 
part of the condensation energy comes from hopping. We 
stress that our results do not change qualitatively if a 
nearest- neighbor hopping is included (the t—t'—J model) 
since symmetry of 4>(k) does not change. 
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For g a this needs some care. Since the G ari (k, ui) has no 
pole at low frequencies we can replace it by G(k, 0). To do 
the frquency sum we need to put back Fb = ujo/ +^0) 
(see Eq. 10) which acts as a frequency cut-off, and take 
the limit — > infinity at the end. 
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